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Abstract 

Spin correlations in metallic and insulating phases of V2O3 and its derivatives 
are investigated using magnetic neutron scattering. Metallic samples have 
incommensurate spin correlations varying little with hole doping. Paramag- 
netic insulating samples have spin correlations only among near neighbors. 
The transition from either of these phases into the low temperature insulat- 
ing antiferromagnetic phase is accompanied by an abrupt change of dynamic 
magnetic short range order. Our results support the idea that the transition 
into the antiferromagnetic insulator is also an orbital ordering transition. 
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Orbital order and disorder have been popular ingredients for explanations of many curious 
phenomena in solids, most notably the transition metal oxides |]],@]. In spite of large interest 
among condensed matter physicists, the evidence for non-trivial, in the sense of not simply 
being due to conventional spin-orbit coupling, orbital fluctuations and order is generally quite 
indirect. We have been able to make progress by taking advantage of a very sensitive measure 
of electronic orbitals, namely the exchange interactions linking the spins on neighboring 
atoms, which are determined by integrals over outer electron wave functions. In particular, 
we report the discovery of a dramatic switching of magnetic exchange interactions on moving 
between the high and low temperature insulating phases of the classic transition metal 
oxide V2O3. The switching is most naturally explained in terms of an orbital ordering 
transition. In addition, we find that the high temperature spin fluctuations in metallic 
V2O3 are insensitive to hole doping, showing little evolution on going from a stoichiometric 
sample with an insulating antiferromagnetic ground state to a sample which remains metallic 
down to the lowest temperatures and displays a low moment spin density wave (SDW). 
Furthermore, the magnetic fluctuations in the paramagnetic metal are similar to those in 
the paramagnetic insulator, which suggest that orbital fluctuations play an important role 
in the metal as well as in the insulator. 

Fig. 1 shows the phase diagram of V2O3 as a function of temperature and two kinds 
of doping. Vanadium deficiency drives the material metallic by populating the 3d band 
with holes, while Cr 3+ (3d 3 ) substitution for V 3+ (3d 2 ) stabilizes an insulating state ||. 
Three single crystal samples studied in this work, and indicated by color bars on the figure, 
cover all four phases of V2O3. The samples were grown using a skull melter and in the 
stoichiometric and vanadium deficient samples, oxygen content was controlled to within 
Sy=0.003 by annealing sliced crystals for two weeks at 1400°C in a suitably chosen CO-CO2 
atmosphere J7|. The sensitivity of our experiment was increased by mutually aligning several 
crystals so that the total mass reached ~15g for V1.973O3, 8g for V2O3 (Neel temperature 
T N =170 K) and 6.4g in the case of (Vo.grCro.os^Og (Tjv=180 K). 

The inelastic neutron scattering experiments were carried out on thermal neutron triple 
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axis spectrometers at the NBSR and HFIR of NIST and ORNL respectively. A py- 
rolytic graphite filter was used to remove high-order neutrons. After correcting for the 
fou;- dependent efficiency of the spectrometer, the magnetic neutron scattering intensity was 
normalized to inelastic scattering from transverse acoustic phonons to yield absolute mea- 
surements of the dynamic spin correlation function [|8]], 

5(Q,w) = \ £(5 Q/3 - Q a Qp)\nQ)\ A j^- 

a/3 

/^e^Ee- jQ(R - R,) (^(t)^,(0)). (1) 
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Wave- vector, Q, will be indexed in the hexagonal reciprocal lattice with a* = 47r/y / 3a = 
1.47(1)A~ 1 and c* = 2vr/c = 0.448(1)A -1 in the metallic phases, a* = 1.46(1)A" 1 and 
c* = 0.449(2)A _1 in the AFI phase, and a* = 1.45(1)A" 1 and c* = 0.451(2)A~ 1 in the PI 
phase. 

We begin by surveying the spin fluctuations in three of the phases of V2O3: the metallic 
antiferromagnet (SDW), the metallic paramagnet (PM), and the insulating paramagnet (PI) 
(Fig. 1). Other workers have successfully surveyed spin waves in the fourth phase, the insu- 
lating antiferromagnet (AFI) ||. Fig. 2(a) and (b) show contour maps of the dynamic spin 
correlation function «S(Q,u;) in metallic antiferromagnetic V1.973O3 at 1.4K and in param- 
agnetic, but still metallic V2O3 at 200K. In both cases, we observe ridges, with bandwidths 
exceeding 18meV and centered near Q=(l,0,0.3) and (1,0,2.3), wave vectors which charac- 
terize the magnetic order in the hole-doped material (e.g., V1.973O3) 0. As is appropriate, 
given the higher temperature for frame (b) than frame (a), the ridges are sharper for frame 
(a). On moving to the PI phase (Fig. 2(c)) at nearly the same temperature, however, there 
is further broadening. Indeed, the data are now best described as a single broad ridge along 
the huo-axis centered at Q=(l,0,0.8). From the half-width-at-half-maximum of constant en- 



ergy cuts through this ridge (see Fig. 3(d) and ref. |]TD[), we estimate spin correlation lengths 
£ c w 1.5A and £ a ~ 2.0A along the c-axis and in the basal plane respectively. 

One of the most remarkable features of the metal-insulator transition in V2O3 is that the 
antiferromagnetic order (see inset of Fig. 1) which develops in the insulator is different from 



the SDW which occurs in vanadium deficient samples |]6|JTT[]. Specifically, as may also be seen 



in the inset of Fig. 1, vanadium atoms in V2O3 have three nearest neighbors within a puckered 
honeycomb plane. In the SDW phase each spin is approximately antiparallel to all of its three 
in-plane neighbors, whereas in the AFI the three-fold symmetry is broken with one nearest 
neighbor parallel, the two others antiparallel. The two types of local spin arrangements 
are conveniently labeled (l(M)-type and (0.5, 0.5, ^)-type respectively, according to which 
line in reciprocal space contains the magnetic Bragg peaks of the corresponding long range 
ordered structure. In the following we show that irrespective of whether we consider static or 
dynamic properties, (0.5, 0.5, £)-type correlations exist only in the AFI phase while (1, 0, £)- 
type correlations exist only outside this phase. 

To probe dynamic correlations corresponding to the two structures, we performed 
constant-energy scans along each of these two directions in reciprocal space. The results 
are shown in Fig. 3 where the right and left columns probe (10£)-type and (0.5, 0.5, £)-type 
correlations respectively. Filled symbols correspond to high temperature and open sym- 
bols to low temperature phases. For both V2_ y 3 (top frames) and V1.g4Cro.06O3 (bottom 
frames), (10£)-type correlations are visible only outside the AFI phase whereas (0.5,0.5,£)- 
type correlations can be seen only in the AFI phase. Entry to the AFI phase not only changes 
the near neighbor correlations, it also brings about coherence in the magnetic excitations as 
evidenced by the resolution-limited double peaks in Fig. 3(a) and (c). These correspond to 
the excitation of counter-propagating spin waves in the long range ordered antiferromagnet 
0. Such dramatic modifications of spin dynamics indicate that exchange interactions un- 
dergo profound changes at the transition to the AFI. At the same time, doping to produce 
either an insulating phase by Cr substitution or a more metallic sample by decreasing the 
V to O ratio has a much smaller effect on the spin dynamics at fixed temperature within 
the paramagnetic phase. Indeed, Fig. 3(b) shows that at 200K, the magnetic fluctuations 
in V2O3 and V1.97O3, which have AFI and metallic SDW ground states respectively, are 
identical. 

Fig. 4 gives the detailed temperature dependence of dynamic spin correlations for V2O3 



and V1.94Cro.06O3- Coincident with the transition to the AFI phase (vertical dashed lines) 
is an abrupt switch between the two types of dynamic spin correlations. A remarkable 
similarity exists between the transitions to the AFI phase from the paramagnetic metal 
(V2O3, left column) and from the paramagnetic insulator (V1.94Cro.06O3, right column), 
which suggests a common mechanism which is independent of whether or not a Fermi surface, 
associated with metallic behavior, exists in the high temperature phase. 

Our discoveries find no comprehensive explanation within a one-band Hubbard model 



l4[13|l . The most serious difficulty is the abrupt switch of the magnetic wave vector which 
occurs at the transition to the AFI. In addition, the insulating paramagnet is characterized 
by magnetic correlations which are shorter ranged than those of the paramagnetic metal. 
Like many single-band Hubbard calculations, the doped copper oxides which eventually 
become high temperature superconductors display none of these anomalous features. 

The most obvious potential cause for the anomalous short range correlations in the PI 
phase is that conventional exchange interactions between spins in insulating (V .97Cr .03)2O3 
lead to a frustrated Heisenberg model which is unable to develop long range order at T~J/kg. 
Solid state chemistry, however, speaks against this possibility because only A and B types 
of nearest neighbors (see Fig. 1) with direct cation-cation overlap have appreciable exchange 
interactions in V2O3 || and a Heisenberg model with only these interactions is not frustrated 



A more likely explanation of our data is based on work of Castellani et al. who estab- 
lished, almost twenty years ago, how covalent bonding between doubly degenerate 3d orbitals 
control the sign and magnitude of spin-spin interactions in V 2 3 0. The basic idea is that 
the magnetically active electron on a single vanadium ion can exist in one of two degenerate 
orbitals. Even so, the exchange integrals, which determine the coupling between spins on 
neighboring ions, depend strongly on which orbitals are occupied. Indeed, when orbitals on 
neighboring ions are orthogonal, the resulting spin coupling is ferromagnetic, while if they 
are not, the coupling tends to be antiferromagnetic. The net Hamiltonian for insulating 
V2O3 then involves orbital as well as spin degrees of freedom at each site, where the spin- 



orbit coupling is not a bilinear coupling on a single site, as in conventional magnets, but 
actually involves the relative spin and orbital occupancies on neighboring sites |2j . Castellani 
and coworkers examined many possible orbital configurations for V2O3 and concluded that 
the most likely is that represented by the colors used to locate the V ions in Fig. 1. Here, 
the magnetic unit cell is doubled in the basal planes because the orbital ordering doubles the 
unit cell also. On warming through TV, we are left with disorder in both spin and orbital 
occupancies. The magnetic short range order in the paramagnetic phase then depends on 
whether we are looking at frequencies (a) large or (b) small compared to the relaxation rate, 
T or b, characterizing the orbital fluctuations. When (a) uj ^> r or b, the interactions with the 
neighbors will be similar to those below T^v, and we expect to see heavily damped renditions 
of the spin waves seen below Tjv- On the other hand, if (b) u <C T orb , one can average over 
the orbital degrees of freedom and obtain an effective spin Hamiltonian where all couplings 
are antiferromagnetic. 

We now consider whether limit (a) or (b) is more appropriate for our (V .97Cr .03)2O3 
sample. Obviously, because we see no remnants of the low-temperature spin waves, descrip- 
tion (a) cannot apply. We can check whether (b) applies, especially in the detailed sense of 
its prediction that all interactions are antiferromagnetic, by rewriting the observed structure 
factor <S(Q, oj) as the Fourier transform of the two-spin correlation function in real space and 
then checking the signs of the near- neighbor correlations. A remarkably good fit (x 2 = 1.5) 



to the experimental <S(Q,u;) [15] is obtained when we truncate the series to include correla- 



tions among only four types of spin pairs (see inset of Fig. 1). More specifically, we find that 
(S • S [0 ,o,i/ 6+5 ]) A = 0.6(3), (S • S [1/3i2/3>I] ) B = -0.19(8), (S • S [2/3 ,i A5 -i/ 6 ]) C = 0.18(8) and 
(S • Sp/3,1/3,1/6]) 13 = —0.09(3), where 5 = 0.026 and we normalized so that (S ■ S ) = 1. 
The curve through the 205K data in Fig. 3(d) was calculated using this model. Thus, 
measurable correlations are not all antiferromagnetic and we do not have a state (b) with 
complete orbital disorder, i.e. with T or b/uj — > 00. Simultaneously, though, the antiferro- 
magnetic correlations are very short-ranged in spite of the fact that the temperature and 
measuring frequencies are below the exchange constants characterizing the AFI phase 0. 
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The most probable cause is that while fcgT and Hu are somewhat lower than T or b, they are 
close enough to r orb that the spin couplings, fluctuating with the orbital occupancies, have 
effectively random signs which give rise to magnetic frustration. In other words, we are 
proposing that the orbital fluctuation rate is of order /cBT^20meV. 

Finally, what happens in the metallic state? Here we can no longer speak about local 
orbital and spin degrees of freedom, but it is still possible to make Wannier projections 
of analogous objects from the band states. Our finding that the magnetic correlations in 
the metallic state are much more similar to those in the paramagnetic insulator than in 



the antiferromagnetic insulator corroborates the assignment by Takigawa et al [16[ of the 
nuclear magnetic relaxation primarily to orbital rather than spin fluctuations in metallic 
V2O3. Even so, the orbital fluctuations may be somewhat faster than in the PI because the 
antiferromagnetic correlations are further ranged (i.e. there are resolvably sharper peaks in 
S(Q,u) for the PM phase) and the PM-AFI transition appears more strongly first-order 
than the PI-AFI transition . 

To summarize, we have discovered that the magnetic fluctuations in V2O3 switch dra- 
matically as a function of temperature as the antiferromagnetic insulating state is entered 
from either insulating or metallic paramagnetic phases. They change in a much more modest 
fashion as a function of doping and temperature in the paramagnetic phase, be it metallic 
or insulating. We conclude that the primary order parameter for the AFI phase is orbital, 
and that orbital order drives the spin ordering not via a conventional spin-orbit interaction 
for each V ion, but instead via a long-range ordered modulation of the exchange constants 
coupling spins on different V ions. The orbital ordering breaks translational symmetry, 
and is the orbital analog of antiferromagnetism. This is to be contrasted with the cele- 



brated orbital order in the cubic uranium pnictides [IT], which while it leads to spectacular 
anisotropies, does not break translation symmetry and indeed can be traced to single-ion 
spin-orbit coupling. Finally, our data represent strong evidence for the long-standing notion 
that orbital, charge, and spin degrees of freedom need to be considered on an equal footing 
near the metal-insulator transition in generic transition metal oxides. 
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FIGURES 

FIG. 1. (color) Phase diagram of V2- y 03-(Vi_ x Cr x )203 in the temperature-composition plane 
[HQ. PI denotes the paramagnetic insulating phase, PM, the paramagnetic metallic phase, AFI, the 
antiferromagnetic insulating phase, and SDW, the metallic spin-density-wave phase. The PM-PI 
phase boundary terminates at a critical point. The color bars mark the samples and the tempera- 
ture ranges explored in this work. Insert: spin and orbital orders in the AFI phase [^Q. The two 
degenerate (in the single ion limit) vanadium orbitals are represented by red and blue respectively. 
A-D label four kinds of near neighbor pairs. 

FIG. 2. (color) Contour maps of S(Q,u) for Q along the (10£) direction in (a) V1.973O3 at 
T=1.4K, (b) V 2 3 at T=200K, and (c) (V .97Cro.o3)20 3 at T=205K. Intensity is indicated by the 
color bars in units of /Mg/meV per unit cell, which contains 6 formulas units. The Q range covers 
a Brillouin zone, with nuclear Bragg points (102) and (104) at the ends. Resolutions in both Q 
and u are narrower than the widths of features by at least a factor 2 [10,llJ (see Fig. 3(b) and (d) 
for a few examples). 

FIG. 3. (color) Constant fiu) = 9 meV scans along (10£) and constant hu> = 12 meV scans along 
(0.5,0.5,£) for V2- y 03 ((a)-(b)) and (Vo.97Cro.o3)203 ((c)-(d)) inside (open symbols) and outside 
(filled symbols) the AFI phase. (10£) scans look similar over explored Tlw-range (See Fig. 2). 
The horizontal bars indicate the FWHM of the projection of the resolution function on the scan 
direction. Peaks in (a) and (c) are resolution limited. We use Ef = 13.7 meV. Collimations are 
60'-40'-40'-80', 60'-40'-80'-80' and 60'-40'-40'-60' at NIST for (a), (b) and (d) respectively; (c) 
shows data from HFIR with 50'-40'-40'-70'. 
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FIG. 4. (color) Temperature variations of the neutron scattering intensity from (lO^)-type 
and (0.5, 0.5,£)-type magnetic fluctuations in pure V2O3 ((a)-(b)) and (Vo.97Cro.o3)203 ((c)-(d)). 
(10£)-type fluctuations at Tiu = 9 meV were probed with Q=(l,0,2) in (a) and Q=(l,0,0.8) in (c). 
(0.5, 0.5, £)-type fluctuations at hu = 12 meV were monitored at Q=(0.5,0.5,0.2) in (b), while we 
report the (0.5,0.5/) Q-integrated intensity in (d). The dashed lines indicate the location of the 
phase transition out of the AFI as determined by the disappearance of a magnetic Bragg peak at 
Q=(0.5,0.5,0). All data were acquired upon increasing temperature from T=1.4K. 
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